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Abstract 

Based on molecular optics we investigate the reflection and refraction of an electromagnetic 
wave between two semi-infinite anisotropic magnetoelectric materials. In terms of Hertz vectors 
and the principle of superposition, we generalize the extinction theorem and derive the propagation 
characteristics of wave. Using these results we can easily explain the physical origin of Brewster 
effect. Our results extend the extinction theorem to the propagation of wave between two arbitrary 
anisotropic materials and the methods used can be applied to other problems of wave propagation 
in materials, such as scattering of light. 
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I. INTRODUCTION 



In classical electromagnetism there are two well-known approaches to the propagation 
of electromagnetic waves. The first is to solve Maxwell's equations with boundary condi- 
iions and the second is to use Ewald-Oseen extinction theorem of molecular optics theory 
^. Compared to the former used traditionally, the latter can give much deeper physical 
insights into the interaction of electromagnetic wave with material Using 



the extinction theorem, the propagations of electromagnetic waves through a semi-infinite 



isotropic material P, 1^ and an isotropic slab [7] have been studied. Recently the Brewster 

mechanism is explained for light incident onto an isotropic material with negative index . 

The extinction theorem also plays a key role in light scattering theory ^. 

In the previous works using the extinction theorem, most deal with the propagation 

of electromagnetic waves incident from free space into isotropic materials. However, the 

opposite situation from an material into free space and especially the situation between 

two materials is rarely studied. Then, whether the methods used in the previous work are 

applicable to the two situations? To answer this question, it is necessary to generalize the 

extinction theorem to investigate the propagation between two materials. On the other 

hand, a recent advent of artificial materials, named as n egat ive-refraction materials, arouses 

^ ^ 

many interest in the field of optics [10|, 
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Since the negative-refraction materials are actually anisotropic magnetoelectric, it is also 
necessary to extend the molecular optics from isotropic materials to anisotropic materials. 

It is the purpose of this letter to use molecular optics to investigate systematically the 
propagation of electromagnetic wave between two semi-infinite anisotropic magnetoelectric 
materials. In terms of Hertz vectors and the principle of superposition we derive the proper- 
ties of propagation and generalize the extinction theorem, so that the propagation between 
two arbitrary materials can be investigated in a unified framework of molecular optics. We 
also explain the mechanism of Brewster effect. The results extend the conclusions about 
the propagation ,n .sot.opic .natenals BBQ and th.ow on tKose obtained by using 
Maxwell equations Q|. The methods used do not require boundary conditions, but can 
reveal the interaction of light with materials, and avoid the difficulties in integrations and 
complex calculations usually encountered in using extinction theorem P|. So the methods 
here can be applied to other problems of wave propagation in materials, such as scattering 
of light. 



2 




FIG. 1: Schematic diagram for the reflection and refraction of TE waves at the interface between 
vacuum and an anisotropic material. The dashed line denotes the possible ray of negative refraction. 

II. REFLECTION, REFRACTION, AND EXTINCTION THEOREM 

In this section, we first employ the formulation of Hertz vector and the principle of 
superposition to deduce the radiated fields generated by dipoles in the propagation of waves 
between two anisotropic dielectric-magnetic materials. Then we derive the real reflected 
and transmitted fields and the Fresnel's coefficients. At the same time, we generalize the 
Ewald-Oseen extinction theorem. Throughout the paper SI units are used. 

Let us consider a monochromatic electromagnetic field of Ej = E^o exp (ikj ■ r — iut) and 
Hj = Hjo exp (ikj ■ r — iut) incident from an anisotropic material into another one filling the 
semi-infinite space z > with kj = ki^ic — ki^z. The schematic diagram is in Fig. ^ Since 
the material responds linearly, all the fields have the same dependence of exp {—iut) which 
will be omitted subsequently. For simplicity, the permittivity and permeability tensors of 
materials are assumed diagonal simultaneously in the principal coordinate system, Sj = 
diag[ej^,ejy,ej,], fij = dmg[fij^, njy, fij,], j = 1,2. 

Following the molecular optics theory, a bulk material can be regarded as a collection of 
molecules (or atoms) embedded in the vacuum. Driven by the external field, the molecules 
are brought into oscillations and then secondary waves are generated by the induced dipoles. 
The radiated electric field by dipoles is decided by jl| 

Erad = V(V ■ He) - Sq^Iq-^ - /XqV X — (l) 
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and the generated magnetic field is 



(9^n 911 
Urad = V(V ■ n^) - + eoV X (2) 

Here lie and 11^ are the Hertz vectors, 

ne(r)= /'^G(r-r')dr', (3) 

Um{r) = J M(r')G(r - r')dr'. (4) 

The dipole moment density of electric dipoles P and that of magnetic dipoles M are related 
to the associated field as P = eoXe ■ E, M = Xm ■ H, where the electric susceptibility 
Xe = (^/^o) ~ 1 and the magnetic susceptibility Xm = (M/yUo) ~ 1- The Green function is 
G{r — r') = exp (zfco|r — (47r|r — r'|), where ko is the wave number in vacuum. In the 
first medium the dipoles produce forward waves as well as backward waves. So, we assume 
the associated dipole moment densities have the following forms 

Pi = PiFexp{ikiF -r) + PiBexp{-ikiB -r), 

Ml = MiF exp {ikiF ■ r) + Mib exp {-ikiB ■ r), (5) 

while P2 = P2F exp {ik2F ■ i") and M2 = M2_f exp {ik2F ■ r) in the second material, where F 
and B label the forward and backward propagating waves, respectively. We now determine 
the Hertz vectors. We firstly represent the Green function in the Fourier form j^. Then, 
inserting it into Eq. Q and using the delta function definition and contour integration 
method, the Hertz vectors can be evaluated 
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(7) 



qjF = kjFx^ + QjzZ, cijB = kjFx^ - QjzZ, q% = k^ - 

biF = kiBx^ + hzZ, hiB = kiBx^ - biz2, = k^ - klBxJ = 1, 2. 



(8) 
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Following the superposition theory of field, the fields in the right region produced by the 
dipoles in the two media add up to the transmitted field. Then, we have 

Et = + El,, (9) 

where the contribution from the first medium to the right side, E^^*^^^*, can be calculated 
by substituting Eq. (jH)) into Eq. ([T)) 

j^l.right _ Q(qiF, Pif) Q( — bjg, Pib) . s 

the field radiated by the second medium itself, E^,, can be obtained after inserting Eq. ^ 
into Eq. (P) 



^2 ^ ^ Q(k2F,f2F) ^ Q(q2F,f2F) ^^^^ 

™ ~ 5'2F 2g2z(92z — k2Fz) ' 



and Q is an auxiliary function 

Q(K, P) = -— [K X K X P + {K^ - eofiou;^)P + cu/xo^oK x M] exp (zK ■ r) (12) 

where M = ■ {k x [P/{eoXe)]/fJ'}/io. Substituting Eqs. ^ and (HD) into ©, and 
considering the transmitted field with the form 

p 

Et = — — exp (iksF ■ r), (13) 

we come to the following conclusions by a self-consistent analysis, (i). From terms with the 
phase factor exp {ik2F ■ i") in Eq- © yield k2F = and the dispersion relation 

1.2 k2 1.2 h.2 



i^2z^2y i^2x^2y ^2zf^2y ^2xP-2y 

for TE and TM waves, respectively, (ii). We know that qi^, q2F and bi^ are all vacuum 
wave vectors. Since only ]i.2F appears in the final transmitted field, they all should be 
extinguished. So we conclude that qiz = q2z = and kip^ = k2Fx = —kiBx- Then, comes 
naturally the Snell's law: kipsinOi = fc2Fsin6't. At the same time, 

Q(qiF, P2f) _ Q(qiF,fiF) _ Q(qiF,fii?) ^ Q (Ig) 

Qlz — k2Fz Qlz — kiFz Qlz + kiFz 

Equation (fT3j) is the generalized expression of the extinction theorem about the forward 
vacuum waves. It describes how the radiation field produced by the dipoles of the first 
medium is extinguished by the counterpart in the second medium. 



Now we study the fields in the first materiaL Similarly, all the fields radiated by the 
whole space are superposed to form the incident and reflected field. 

E, + E, = El, + Eltf, (16) 

where the field radiated by the first medium, E^^^, can be calculated substituting Eq. © 
into Eq. jH) 

the contribution E^^'^-^* from the second medium to the left half-space can be obtained after 
inserting Eq. ((7j) into Eq. ((H) 

T^2.left ^ Q(q2B, P2f) 

™' 2q2z{q2z + k2Fzy ^ ^ 

The incident and reflected fields may be written as 

Ej = exp {ikiF ■r),Er = exp (-ikiB ■ r), (19) 

respectively. Inserting Eqs. (fTTj) . (fTK|) and (fTIHl into which must hold true for everywhere 
in the left half-space, we come to the following conclusions, (i). From terms with the phase 
factor exp {ikip ■ r) or exp {ikis ■ r) in Eq. (fTBj) follow that kip = kj, ki^ = k^. and the 
dispersion relation is like Eq. ()14|1 after replacing the subscripts 2 with 1 and t with i, 
respectively, (ii). qi^ and q2B are vacuum wave vectors and should be extinguished. So, 
we have qie = q2_B and 

Q(qiB,PiF) ^ Q(qii?, Pib) _ Q{q.iB,P2F) ^ Q ^20) 

qiz + kiFz qiz — kiFz qiz + ^2Fz 

which is a new expression of the extinction theorem we find. It shows how the the backward 
vacuum field produced by the dipoles of the second medium is extinguished by that in the 
first medium. 

Solving the set of Eqs. (fT^. (fT3j) . (fT^. and (j^ . we obtain the reflection coefficient 
Re = Ero/EiQ and the transmission coefficient Tg = Eto/EiQ for TE waves 

^E = -J : J^E = -J : {^i-) 

fJ'2xKiz + jJ^lxl^tz l^2xi^iz + l^lxl^tz 

Analogously, we can discuss TM waves and obtain similar results after replacing /x in with 
£, considering the expressions of E in Eq. ((H) and H in Eq. ((2)). 
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Hence, we have obtained the generahzed extinction theorem, i.e. Eqs. fll5|) and ()20p. 
So, the propagation between two arbitrary materials can be studied in a unified framework: 
Under the action of external field in the first material, the molecules in the tow materials 
are driven to oscillate and generate induced fields. The transmitted wave is the result of 
superposition of the incident vacuum field from the first medium and the fields radiated by 
the induced dipoles in the second medium. The refiected wave is the sum of the backward 
vacuum field from the second medium and the backward fields induced by the dipoles in the 
first medium. Note that Eqs. (fT^ and (plj) should hold true for every point in the associated 
half-spaces, which indicates that the cancellation of vacuum waves occurs everywhere inside 
the materials. 



III. ORIGIN OF BREWSTER EFFECT 



In what follows we apply the above conclusions to discuss the origin of Brewster effect. 
If the power refiecti vity \R \^ = 0, there is no refiected wave and the incident angle is called 
Brewster angle [s^l^QjHl- From Eqs. ()20|) follows the refiected field magnitude 

1 fqiB X [qiB X {xie ■ Eio)] , qiB X {xim ■ [/ioA^r^ ■ (ki X Eio)]} 



qiz + kizHo/fJ'ix { qiz + kiz ' qiz + h 

^ qiJ? X [qijj X {X2e ■ Efo)] ^ qij; X {X2m " [/ioAt^^ ■ (kt X Effl)]} \ 

qiz + ktz qiz + ktz J ' 

On the right-hand side of Eq. (j221), the first two denote the contributions (labelled as E*q) 
of dipoles in the first medium to the refiected field, while the last two are the contributions 
(labelled as E*q) from the second medium. In order for zero refiection, it requires that 
E*Q + E*Q = 0, from which follows the condition for Brewster effect: If 



fJ'2z£2y — f^lz^ly 

the Brewster angle for TE waves is 



> 0, (23) 



qTE ^ ^^^-1 f^Uf^2ze2y-f^izeiy) _ ^24) 

y f^lzf^2z{f^2x£ly — fJ'lx£2y) 
If /i2a7/ilx = £2y/£ly^£2y/£ly 7^ fJ'lz/fJ'2z, then 9]^^ = 0] When /i2z//ila' 7^ £2y/£ly^£2y/£ly = 

fi'iz/fJ'2z, &W = ^/2; If fJ'2x/fi'ix = £2y/£iy ^ ^2yl<^\y = fJ'iz/fJ'2z, then Brewster effect will 
occur for any angle of incidence, which may lead to important applications in optics. All the 
typical and nontrivial sign combinations of e and fj, are shown in Table. lU Note that, if the 



TABLE I: Existence conditions of Brewster angles for TE waves at the interface of two anisotropic 



media. 
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Note: * denotes either of it, V indicates that Brewster angle exists for any parameter values. 

first medium is vacuum, E*q = 0, then zero-reflection can occur only when E*q = 0. In other 
words, the Brewster effect is because the contributions from electric and magnetic dipoles 
of the medium to the reflected fleld in vacuum add up to zero, which is also pointed out in 
Ref. 0]. In order to describe these conclusions, we give numerical examples in Fig.|21 We also 
follow the method in Ref. 17| to simulate in Fig. El a beam 

incident from an anisotropic material to another where f{k±) is the Gaussian modulation. 
We can see numerical simulations are in agreement with the theoretical conclusions. Ex- 



perimentally Brewster effect has been realized with metamaterials [26[. So, one can realize 
zero reflection through choosing appropriate material parameters, which may lead to make 
polarizers or light splitters. 
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FIG. 2: Reflectivity and reflected field magnitudes, normalized by the corresponding incident field 
magnitudes, for TE wave incident incident on the interface between two anisotropic materials. The 
three curves of refiectivity in (a) correspond the case in (b), (c), and (d), respectively. Since the 
first medium is vacuum in (b), the radiated electric fields (E*q) generated by the second medium 
forms the refiected field, then Brewster angle 0™ appears when E*q = 0. In (c), 6c is the critical 
angle of incidence. In (d), E*q + E*q = 0, so omnidirectional total transmission occurs. 

IV. CONCLUSION 

In summary, we carried out a systematical investigation on the propagation of wave be- 
tween two anisotropic magnetoelectric materials. In terms of Hertz vectors and the principle 
of superposition we derive the properties of propagation and generalize the extinction the- 
orem, so that the propagation between two arbitrary materials can be investigated in a 
unified framework. We apply the results to explain the physical origin of Brewster effect. 
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FIG. 3: Reflection and refraction of a Gaussian beam incident on the interface between two 
anisotropic materials. For (a)-(c) and (d)-(f) the two materials correspond to those in (c) and 
(d) of Fig. [2 respectively. The phase and the energy flow are refracted regularly and anomalously, 
respectively in (a)-(c), but are both refracted anomalously in (d)-(f). 

The methods do not require boundary conditions, but can disclose the microscopic process 
of light propagation in materials, and avoid complex calculations usually encountered in 

wave 



using the extinction theorem. So the methods can be applied to other problems o 
propagation, such as scattering of light ^j, propagation through stratified media 3| 



metamaterials l29|, and so on 



or 
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